Branching of negative streamers in free flight 
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We recently have shown that a negative streamer in a sufBciently high homogeneous field can 
branch spontaneously due to a Laplacian instability, rather than approach a stationary mode of 
propagation with fixed radius. In our previous simulations, the streamer started from a wide initial 
ionization seed on the cathode. We here demonstrate in improved simulations that a streamer 
emerging from a single electron branches in the same way. In fact, though the evolving streamer is 
much more narrow, it branches after an even shorter propagation distance. 



Streamers are widely known phenomena in electric 
breakdown: long and rapidly growing channels of high 
ionization penetrating a region of low or vanishing ioniza- 
tion under the influence of a strong electric field 0,^. In 
experiments, streamers frequently are seen to branch. On 
a phenomenological level, the so-called dielectric break- 
down model (DBM) |^-|| , a variation of diffusion limited 
aggregation has been suggested as a stochastic pro- 
cess explaining the branching. Such models are inspired 
by earlier concepts of streamer propagation where the 
space charge was assumed to be smeared out over the 
full streamer head. Branching then would occur due to 
randomly distributed ionization avalanches around the 
streamer. Such stochastic phenomenological models on 
larger scale have not been related to particular discharge 
models. 

However, simulations |^,^ have established that the 
space charge of the streamer is concentrated in a thin 
layer around the head rather than being smeared out over 
the full head. This leads to a different field distribution 
and a much faster propagation mode, and to a revival of 
the concept of an ideally conducting streamer formulated 
by Lozansky and Firsov In particular, this charge dis- 
tribution allows for a completely different "Laplacian" in- 
stability mechanism of streamers 0Jll[ that can operate 
even without any randomness quite like in viscous finger- 
ing, dendritic growth of solids into undercooled melts etc. 
It naturally emerges in the minimal streamer model (|l|) 
~ (||) . In a recent paper jlj] , we found that in contrast to 
previous expectations, this Laplacian instability mecha- 
nism can cause the spontaneous branching of streamers 
propagating in a strong homogeneous field. The system 
was identical to that of previous simulations |^J^ except 
that the field was twice as high. However, our simulation 
results were somewhat limited by numerical constraints. 
Because of the unexpected results and some non- 
smooth structures in the figures, some researchers won- 
dered whether the figures showed a numerical rather than 
a physical instability Q . 

Therefore we here present new simulations with im- 
proved numerics, and we show more details of the evo- 
lution. In particular, we increased the number of grid 
points in the simulations from 1000x1000 to 2000x2000, 
which improves both the accuracy of the numerical out- 



put data and the quality of the contour plots, so the new 
plots are smooth. Furthermore, we changed the bound- 
ary condition on the cathode from homogeneous Neu- 
mann d^a = to homogeneous Dirichlet a — for the 
electron density a. This means that while previously the 
electrons freely could fiow from the metal of the electrode 
into the gas, this current is now suppressed. With these 
improvements, the choice of the initial conditions now 
ceases to be constrained by numerical considerations. We 
continue to use cylinder symmetry to calculate effectively 
on a 2-dimensional grid. Doing so, we keep relying on 
the analytical argument that the constraint of cylindri- 
cal geometry suppresses some instability modes, and that 
therefore a truly 3-dimensional system would become un- 
stable earlier or at the same time as the system with 
symmetry constraint, but certainly not later. 

The improvements of the numerical code allow now for 
the simulation of a streamer starting from a single elec- 
tron on the cathode rather than from a wide initial seed 
as in I p^ . Our new figures show the streamer branch 
in free flight. The initial electron on the cathode cre- 
ates first an avalanche, then a streamer, and finally the 
streamer splits. The situation resembles the historical ex- 
periments of Raether |jl[, except that the field is higher. 

In detail, we investigate the minimal streamer model, 
i.e., a "fiuid approximation" with local field-dependent 
impact ionization reaction in Townsend approximation 
Q in a non-attaching and non-ionized gas. In dimen- 
sionless units 00-^3], the model has the form: 



dta - V • (fj E 



D Vcr) 


-^/(|E|) , 


dt p 


-^/(|E|) , 


p-a 


= V • E , E = -V$ 


/(|E|) 


= |E| e-i/|E| . 



(1) 

(2) 
(3) 
(4) 



Here a and p are the densities of electrons and posi- 
tive ions, E is the electric field, $ the electric poten- 
tial and D the dimensionless diffusion constant. The 
mobility of ions is neglected. The same model was in- 
vestigated in @J^. The translation to physical units 
depends on the type and the pressure p of the gas. 
With an effective field-dependent impact ionization coef- 
ficient a{E) = Ap exp (— i?p/|£'|) and electron mobility 
/ie = ft/p as in P], and with the parameter values for 
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nitrogen as in [0,D, we get for the scales of length, time, 
field and particle density 
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where [N2] is the neutral gas particle density. The pres- 
sure dependent scaling relations are of interest for labora- 
tory experiments, as well as for so-called sprite discharges 
||15| in the mesosphere above thunderclouds. 

The simulation is performed with the same voltage of 
1000 and electrode distance of 2000 as in jlj], correspond- 
ing to a field of 100 kV/cm and an electrode distance of 
4.6 mm for nitrogen at p = 1 bar. As an initial condition, 
one electron was distributed over one cell of the numer- 
ical grid next to the cathode. Because of the diffusive 
broadening of the drifting avalanche, this "distribution" 
of the initial electron is permissible. 

In the figures, we omit the avalanche phase and con- 
centrate on the developed streamer and its branching at 
times t = 500, 525, 550 and 575. The field is apphed 
in the z direction, and r is the radial coordinate extend- 
ing up to r = 2000 to suppress artifacts from the lateral 
boundaries. The simulation shows a streamer not con- 
nected to the electrode. It has the conical shape of the 
initial avalanche created by a single electron in a homo- 
geneous field as first found in Raether's experiments ||l|]. 
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FIG. 1. Evolution of an anode directed streamer in a strong 
homogeneous background field. The planar cathode is located 
at z = and the planar anode at z = 2000. The radial co- 
ordinate extends from the origin up to r = 2000 to assure 
homogeneous field conditions. Shown is 100 < z < 750 and 
< r < 150 with equal axis scaling. The thin lines denote 
levels of equal electron density a with increments of 0.15. 

Fig. 1 shows the electron density a within the whole 
streamer body at four snapshots of the evolution. Due to 



the no-current-boundary condition on the electrodes, the 
streamer is not connected to the cathode at z = 0. The 
equidensity lines mark multiples of 0.15 in all snapshots. 
The electrodes at z = and z — 2000 are outside the 
region shown. The qualitative evolution is the same as 
in our previous paper [^ : while the streamer extends, it 
becomes wider, the electron density increases, and finally 
the instability develops. 
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FIG. 2. The streamer at time t — 525. The figures show 
contour plots with equal aspect ratio for the electron density 
a, ion density p, space charge density p — a and electric poten- 
tial <I>. The level lines for the densities a and p are multiples 
of 0.1, for the space charge density p — a oi 0.05, and the 
increment between different levels of the potential is 10. 

Fig. 2 shows four different observables of the whole 
streamer at one time instant t = 525. We plot the elec- 
tron density cr, the ion density p, the space charge density 
p — (T, and the electric potential <I>. In the propagating 
front at the tip of the streamer, an overshoot of drifting 
electrons creates a negatively charged layer, while the 
back of the streamer is depleted from electrons and pos- 
itively charged by the essentially immobile ions left be- 
hind. Therefore the "ion streamer" is longer in the back, 
and the ion density has no overshoot in the propagating 
tip. The resulting space charge densities are shown in 
the third plot: a positive space charge with a maximum 
density of 0.25 in the back, and a negative space charge 
with a minimum density of —0.25 in the propagating tip. 
While the whole object is electrically neutral, the elec- 
tric polarization leads to a suppression of the field in the 
interior. This can be read from the large distance be- 
tween equipotential lines inside the ionized body visible 
in the last plot. The streamer is indeed approaching the 
Lozansky-Firsov limit of "ideal conductivity" ^ . 

Fig. 3 shows a zoom into the tip of the streamer at the 
same time steps as in Fig. 1. Again the electron densities 
a are plotted, and addionally the equipotential lines $. 
Already at time t = 525, the onset of the instability can 
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be concluded from the somewhat irregular shape of the field ahead of the streamer increases. Consequently, the 
level lines. The curvature of the tip decreases, the den- tip becomes unstable and splits, 
sities increase, the field in the interior decreases and the 
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FIG. 3. Zoom into the evolution of the streamer tip showing electron density a with increments of 0.1 as thin lines and 
electric potential $ with increments of 5 as thick lines. At the last time step t = 575, the maximal electron density in the core 
of the finger on the axis is 1.5. AH figures have equal aspect ratio and equal axis scaling. 



We conclude with a few remarks: 

1. It is quite remarkable that the finger of the new sim- 
ulations branches after an even shorter travel distance 
than the one growing out of a wide, highly ionized seed in 
[^2|, though it is more narrow during its whole evolution. 

2. When the initial electron is not placed on the electrode 
but somewhere in the gap, simulations show the same 
figures as in this paper, essentially unchanged except for 
a trivial shift in space. 

3. The transition times from avalanche to streamer and 
then further from streamer into the instability depend 



quite nonlinearly on the externally applied field as we 
will discuss in a forthcoming paper; therefore the ques- 
tion whether branching also would occur in a lower field 
after a longer propagation distance, is open. 
4. Finally, the discovery of sprite discharges in the 
mesosphere gives a new impetus to the study of negative 
streamers in high fields, since a negative sprite discharge 
propagates into increasing heights, i.e., areas of decreas- 
ing pressure; the scaling laws (^) therefore imply that 
the effective field continues to grow along their path. 
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